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1. Introduction 



There are a wealth of apphcations of 1 + 1-dimensional quantum field theories defined 
upon the half-line. Such theories generically have some potential term at the boundary, and 
are thus often referred to as "boundary field theories." Amongst other things, they have 
been used to describe dissipative quantum mechanics fll], Kondo effects 0, and quantum 
wires . The exemplar of such theories is the sine-Gordon model, whose Lagrangian is 
given by: 

dx + Mcos^{(f{x = 0)-(po) . (1.1) 

The massless limit for = 8n was considered in [H, and for generic values of (3'^ in [^. 

The quantum version of the sine-Gordon model was considered by Ghoshal and 
Zamolodchikov 0, who conjectured its integrability and an exact S-matrix. Our purpose 
here is to discuss these two aspects in the classical and semi-classical limits, generalizing 
the well known bulk analysis to the boundary problem (|l.lj). 

We establish the classical integrability in the next section by using the bulk sine- 
Gordon theory and the method of images. Although such a method might not be expected 
to work, a priori, due to the non-linearity, it turns out that the boundary condition in- 
herited from (|1 . 1[) can be realized, for the boundary refiection of a single soliton, by using 
a three soliton solution of the bulk problem. We compute the classical time delay for the 
most general values of the parameters M and (po. The semi-classical analysis is performed 
in the third section of this paper. We find agreement with the (3^0 limit of the results 
in0. 



(dtif)'^ - (d^iff + g cos (3(p 



2. The classical solutions 

For a conformal field theory, or for any linear partial differential equation, the simplest 
method of dealing with a boundary is to use some version of the method of images fl^ 



One of the purposes of this section is to describe how this technique can also be applied 
to the classical sine-Gordon equation on the semi-infinite interval. Indeed, we will show 
that this method replicates the soliton scattering from a boundary with the most general 
boundary conditions consistent with integrability 0. 
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2.1. The T -functions 

After appropriate rescaling of the field and the couphng constant we have, for the 
classical problem associated with the bulk part of (|1 . 1|) , the sine-Gordon equation 



= -sin((^) . (2.1) 



where fiip = (p. On the interval [0, oo) it was argued in |^ that the most general boundary 
conditions consistent with integrability were: 

5,0U=o = Msin(i((/>-(/>o))|^^o , (2.2) 

for some constants M and (pQ, corresponding to the boundary term in ( |1 . ![ ) . 

On the infinite interval, (—00,00), the classical multi-soliton solution to the sine- 



Gordon equation is well known [^. It is usually expressed as: 



4>{x,t) = 4arg(r) = 4 art an (^^^ ) , (2.3) 



where the r-f unction for the A^-soliton solution is: 

N 



T 



J2 e^^^.-i^^'^^^ exp -J2'2N 



kj J \ kj 
+ 2 ^ ^i^j log ( 



l<i<j<N 



ki -\- kj 



(2.4) 

The parameters kj, aj and ej have the following interpretations. The velocity of the j^^ 
soliton is given by: 

(Note that Vj is positive for a left-moving soliton.) The aj represent the initial positions 
of each of the solitons, and ej = +1 if the j^^ soliton is a kink, while ej = —1 if it is an 
anti-kink. The rapidity, 6, of the soliton is defined by /c = e^, and we have normalized the 
soliton masses to unity (in further discussion, the words "soliton" and "kink" will be used 
synonymously) . 

It is fairly obvious how to get a single soliton solution on [0, 00) with either (p\x=o = 
or dj:(j)\j:=o = 0. One exploits the symmetry of ( pTl| ) under (p —cp and x —* — x, and 
simply takes a two soliton solution on (—00, 00) where one soliton is a mirror image of 
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the other through x = [12|. If one does this with a double-kink solution then it satifies 
the foregoing Dirichlet condition, while the kink-anti-kink solution satisfies the Neumann 
condition. It is also not hard to guess how one can go beyond this solution: For M = oo, 
the boundary condition ( p72|) reduces to (t)\x=o = (po- The only way that this can be 
obtained from a multi-soliton solution on (—00, 00) is to put a third, stationary soliton at 
the origin. 

We therefore consider the three soliton solution with ki = k, k2 = 1/k and k^ = 1. 
That is, we consider: 

1 1 

T = 1 — €V e k"- ^ ' — gQ ^—i e 2k'- ^ ' Fit) 

(2 6) 

where we have introduced the shorthand: 

e = eie2 , eo = £3 , a = ai + a2 , b = as . (2.7) 
The function F{t) is defined by: 

F{t) = eie-^(^'-^)*-"i + e2e^(^'-^)*-"^ . (2.8) 



This solution has = at x = 00, and for /c > 1 it describes a left-moving soliton moving 
from X = 00 with a right- moving "image" starting at x = —00. There is a stationary 
soliton with center located at x = —6. Viewing this as scattering off a boundary at x = 
one can easily see that a is the phase delay of the returned soliton. To make this more 
explicit, observe that r has the following asymptotic behaviour: 

r(x,t) ~ 1 + ,,^-J^l(k' + l)-+ik'-m-a^ ^,_i^oO with f = -g^; 

^ t fc +1 ^2.9) 

r(x,t) ~ 1 + e-2fc[('='+i)"-(^'-i)*]-"2 x,t^oo, with f = . 

The problem now is to first show that the r-function given by ( p.6| ) provides a solution to 
the boundary value problem on [0, 00) defined by ( pTll ) and (p.2|) . Our second purpose is 



to relate the parameters a and b of ( |2.6| ) to the parameters M and ^0 of (p72|), thereby 
obtaining the classical phase delay, a, in terms of M and (pQ. 



2.2. The classical phase-delay 



To summarize the computation briefly, one substitutes ( p.6|) into ( |2.2|) , and obtains 
the constraint: 

[7^e(r)a^Xm(r) - 7^e(r)a^Tm(r)] 

= M [2cos(i(/)o) 7^e(r)Jm(r) - sin(|0o) (7ee(r)2 - Xm(r)2)] |^^^ . 

(2.10) 

One can solve this by brute force substitution for the real and imaginary parts of r, but it 
is somewhat simpler to find constants a, /3, 7 and S such that: 



lx=0 



(2.11) 



a,7^e(r) I^^Q = [a7^e(r) + /3 Jm(r) 
a.Jm(r) I^^Q = [7 7^e(r) + 6 Im{r)]l^^ 

One then finds that ( g?TO|) can be satisfied if and only if one has a — 6, which indeed turns 
out to be true. Using this one arrives at: 



M cos(i(/)o) 
M sin(i(/)o) 

where 



-2 eo 



1 + 



(2.12) 



(2.13) 



It is algebraically very tedious to invert this relationship. One proceeds by eliminating 
e~^, and then solving for a. It is very convenient to introduce a new parametrization of 
M and 4)q: 

H = M cos(^(^o) = 2 cosh(C) cos(?7) ; 

u = M sin(^(/)o) = 2 sinh(C) sin(77) , 

where < ( < 00 and —n < rj < tx. (In the (/U, z^)-plane the curves of constant ^ are 
ellipses, while the curves of constant 77 are hyperbolae whose asymptotes make an angle of 
with the /U-axis.) We then find that the phase delay is given by: 



(2.14) 



a = log<! -e (tanh (f))^(tanh(^)) 



tanh ^ (6* + ir]) tanh \ {0 — irj) 



±1 



(2.15) 



tanhi(^ + C) tanhi(^-C) 

There are several things to note about this formula, 
(i) The ambiguity of the ±1 power comes from solving a quadratic equation, and is a 
direct reflection of the fact that (|2.2| ) is not invariant under (p —>■ (p + 2n (whereas ( |2.1| ) 
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is invariant under this shift) . Equivalently, one can flip between the + and — powers 
by sending — > (/>o + 27r. 
(ii) The argument of the log is always real and positive. The discrete parameter, e = eie2, 
must be chosen to arrange this. Hence: 

e = +1 for - C < ^ < C e = -1 for 1^1 > C • (2.16) 

This means that a kink reflects into kink for — C < < C, and reflects into an anti- 
kink for 1^1 > C. This is consistent with the fact that Dirichlet boundary conditions 
(M = oo) cause a kink to reflect as a kink, whereas Neumann boundary conditions 
(M = 0) cause a kink to reflect as an anti-kink. Note that these two domains of 
parameter space (in which a kink reflects differently) are separated from one another 
by a logarithmic singularity in the classical phase delay. 



(iii) The choice of the power ±1 in (|2.15|) is correlated with the parameter ry and whether 



there is a kink, or anti-kink at the origin. Speciflcally, we have: 

eo = ± si^n(tan(§)) , (2.17) 

where the ± is the same as that of ( p.l5| ). 
(iv) In the M ^ oo, or Dirichlet, limit we see that: 

C ~ log(M) , ry ~ |(/)o and e = -1 , (2.18) 

and the phase delay collapses to: 



a = log<^ (tanh(f))^(tanh(6'))' 



tanhi(^ + i^) tanhi(^-z^) 



±1 



(2.19) 



In this limit (|2.2|) enforces Dirichlet boundary conditions. It is, however, important 
to note that there are two possible distinct boundary values: (/)\x=o = <Po and (/>|x=o = 
(po + 2ti. Since the boundary potential is —Mcos{^{(p — (j)o)), one sees that (j)\x=Q = 4>o 
is stable, while ^|x=o = + Stt is unstable. 

From now on we consider only the stable boundary value that corresponds to the 
positive sign in (|2.19|) . Then one has from (|2.17|) 



eo = - sign{tfm{\(t)\x=Q)) . (2.20) 

It is essential to observe that we have taken (l)x=oo = 0, ah initio. For different 
boundary conditions at x = oo, one should replace by 4>\x=o — 4>\x=oo- This physical 
parameter is deflned mod A-k. An independent derivation of time delay for the Dirichlet 
boundary conditions, along with further details, are provided in the Appendix. 
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2.3. Boundary breather solutions 



To fully understand the semi-classical scattering computation one also needs another 
class of classical solutions, which we call here "boundary breathers" . It is well-known that 
"breathers" represent bound states in the soliton-anti-soliton channel in the bulk sine- 
Gordon theory. In the same way that the classical bulk breathers can be obtained from 
the appropriate solution by analytic continuation of 9 to imaginary axis, one might expect 
that the same procedure would give boundary breathers on the half-line. To see this, we 
set ^ = (0 < t? < ^) in the formula (Al) (see appendix). Next, we impose the following 
conditions for a solution to be boundary breather: 

a) it should be a real function, 

b) it should have finite energy and 

c) the asymptotic value at x = -|-oo must be fixed and equal to 27m (n - integer number). 
The three-soliton (resp. soliton-anti-soliton-soliton) configuration satisfies the first 

condition provided that "i? < — ^ (resp. "i? < ■^). However, the other conditions are sat- 
isfied by the soliton-anti-soliton-soliton configuration only, which one could have foreseen 
from the analogy with the bulk theory. The boundary breather solution has the form (see 
figure 2): 

(pb = 27r- 

2 cot 1? cot cot ^e^+^ ^ cos(t sin i9) + e^^ "^K cot^ | + 1 (2.21) 

— 4 arctan 1 — = = r ^ „ , 

2 cot d cot f VK ™^ ^ cos(t sin ^9) + cot ^ (e^^ ^i^ + cot^ | ) 

where 

K = cot(— - -) cot(— + -). 

So, we have a continuum of classical boundary bound states when < < < ^ 
and for other (po according to the 27r-periodicity. In the quantum theory this continuum 
shrinks into a discrete set of bound states (see next section). Note that in the limit ^ ^ 
the boundary breather ( p.21|) reduces to the ground state, figure 3, and the phase delay has 
singularities at 6* = ±^(/)o- An analogous picture of bound states occurs for the anti-soliton 
scattering with fixed boundary conditions. 
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2.4- General solutions, integrability and BdcMund transformations 

Thus far we have only apphed the method of images to obtain certain special classical 
solutions of the boundary sine-Gordon problem ( p. . 1|) . It is natural to suggest that general 
solutions of ( p. ■ 1| ) can be obtained by similar methods. This in turn would establish the 
classical integrability the boundary problem ( |1 . 1[ ) . It is, in fact, rather straightforward to 
show that both of these conclusions are true, at least for the problem ( |1 . Ij ) with (po = 0. 
The method we will employ can almost certainly be generalized to problems with (po ^ 0, 
and also has the virtue that it can be used to construct the integrable boundary potentials 
for the more general To da models. A related approach has been followed in The 
basic idea is to use the fact that any integrable hierarchy has Backlund transformations: 
that is, solutions can be mapped into one another by non-trivial gauge transformations 
constructed from the afRne Lie algebra action on the corresponding LAX system. For the 
sine-Gordon equation, the requisite Backlund transformation can be cast in the following 
form: 

d^{<j) + ij) = e'^ sinf^^) ; d^{<j) - ^) = e"^ sin f , (2.22) 



where u = x — t, v = x + t, and ( is an arbitrary constant parameter. The point is 
that (f) satisfies the sine-Gordon equation ( pTTI) if and only if i/j does so as well. Suppose 
that '0 is a solution to sine-Gordon on [0, oo) satisfying a Dirichlet boundary condition: 
ip\x=o = 2?7, where ry is a constant. It follows immediately from (|2.22| ) that (p is a solution 
to sine-Gordon satisfying (|2.2|) , where M and cpo are given in terms of C and r] by ( |2.14|) . 



Thus, if one can solve the Dirichlet problem, one can solve the more general problem by a 
Backlund transformation. 

Observe that if i] = 0, or equivalently (po = 0, then the Dirichlet problem can be 
solved trivially by method of images: one gets the solution on the half-line by extending 
it as an odd function on the full line. Thus, the Backlund transformation essentially 
defines the generalized method of images. It is also by no means an accident that the 
parameters entering into the Backlund transformation (77 and Q are precisely the rapidity 



parameters that turn up in the phase delay (|2.15|) . One should also note that the form of 



the integrable boundary potential is given directly by the Backlund transformation. This 
fact should easily generalize to Toda systems . 

Backlund transformations, in general, are invertible transformations on the solution 
space of an integrable hierarchy. The simplest forms of them modify the constants of 
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the motion of solution, and possibily add or subtract a sohton. One can certainly find a 
Backlund transformation that will introduce a stationary soliton into the soliton-soliton 
solution of sine- Gordon. As a result, the general three soliton solution employed above 
can be obtained from the two soliton solution that is appropriate for the "trivial" Dirichlet 
problem with (po = 0. We therefore expect that any solution of the trivial Dirichlet problem 
can be mapped onto the generic problem ( |1.1D , thus establishing the classical integrability. 
Here we shall content ourselves with having explicitly established these results for = 0, 
and having shown that there is a direct link between the form of the potential and the 
structure of the underlying integrable hierarchy. 



3. The semi-classical analysis 

We start the discussion by considering the jS'^ limit of the results in 0. The 
(quantum) boundary S matrix elements are 



S+{9) = P+{e) = cos - (t - l)i9] R{e) 
SZ{e) = P-{9) = COS + (t - l)i9] R{9), 



Stt 



The function R{0) decomposes as 



(3.1) 



and 

s+ = Q+{e) = s+ = Q-{e) = ^ sin [2(t - i)ie] R{e) , (3.2) 

where 6 is the rapidity of an incoming particle. The parameter t is defined by: 



R{e) = Ro{e)Ri{e), (3.3) 



where Rq is a normalization factor ensuring unitarity and crossing symmetry that does not 
depend on the boundary conditions. The dependence upon boundary conditions appears 
in which reads 

Ri{e) = -^MTi,e)a{iQ,e). (3.4) 

cos^ 

Two of the four parameters /c, ^, ry, O are independent, and we have the relations 



1 . o . 2 



cosrycoshG = — cos^, cos^ ?7 + cosh^ G = 1 + j^. (3-5) 



These parameters are related with M and (po in an unknown way. An expression for the 
functions Rq and a involving infinite products of F-functions is given in 0. There are also 
simple integral representations: 



cos[E, — i[t — 1)0] \J-oo ^ 2 cosh(r — Ijx smhx I 

and 

B„W=exp|-r '-^ ^inhi(f-l)x.inh(i-l). ^.j.^^.A 

^ ^ ^\ X sinhf sinh2(t- l)a; J ^ ^ 

One has cr(^, ^) = cr(— ^,6*). The only difference between the scattering of solitons and 
anti-solitons arises therefore from the pre- factor cos[^ ^ i{t — 1)9] in ( |3.1J ). 

For simplicity we consider only the limit in which M — oo. One can identify the 
corresponding values k = and Q = oo easily since, at these values, the topological charge 
Q = 2^ /o°° ^xfdx is conserved and therefore the amplitudes Q± must vanish. One has 
also ?7 = ^ so 

Ri{e) = ^a{^,e). (3.8) 

COS c 

Consider now the quantum boundary S matrix at the leading order in as (3 0. 
The computation is most easily done by using the integral representation given above, 
and evaluating the integrals explicitly by the residue theorem. This provides convergent 
expressions where the /3 ^ limit can be investigated term by term. To get non-trivial 
results one must scale ^ as 

^^^^ 

We find then 

S{e;y) = exp (^j'dvlntaiih'' j , (3.10) 



where 



S{9; 0) being the semi-classical approximation to the bulk soliton-soliton S-matrix ^ (in 
the following we denote it also as S{9)), and k = sign^ (in the following we assume that 
9 > 0). 

Before discussing the relation between formula (|3.9| ) and the classical computations of 
the preceding section, it is useful to comment on the bound states of the quantum theory. 
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Poles of the Rq term are associated with breathers and do not correspond to the boundary 
(new) bound states. The latter correspond rather to poles of the cr(^,6') term which are 
located in the physical strip Im9 G [0,7r/2]. By inspection of the F-product expression 
one finds two families of poles of a: 

l> (3.11) 



where n, I are integers. Let us restrict to ^ > 0. Then the only physical poles are those 
that correspond to + sign in 6'(i) and — sign in 9^'^\ The first pole that enters the physical 
strip (from the bottom) is ^q^q for ^ > ^. The number of poles of cr in the physical strip 
increases monotonically with ^ for ^ small enough, and as t gets large it becomes simply 
of the form ^/tt. These poles are cancelled by the cosine term in P_ and therefore appear 
only in the P+ amplitude. As ^ reaches the value ^ = the poles densely fill the interval 
< Im6' < ^ and a pole at 6 = z7r/2 appears corresponding to the emission of a zero 
momentum soliton. As argued in ^ this corresponds to a change in the ground state of 
the system. For < cpo < ^ the ground state is — > at infinity but for ^ < (po < 3^ 
it is — > ^ at infinity. Therefore the value ^ = ^p- corresponds to (^o = ^ and it is 



the upper physical value for ^: as (po varies arbitrarily, ^ never gets larger than and 

(2) 

the set of poles 9^^ [ never enters the physical strip. Observe that there are bound states 
in the quantum theory when in the classical case the kink sitting in the middle and the 
incoming one are of opposite topological charges in complete agreement with the discussion 
of boundary breathers in section 2. 

To establish the relation between ( |3.9| ) and ( |2.19| ), first recall [H the general relation 
between the quasi-classical scattering phase shift, 5(6*), and the classical phase shift, a{6) : 

5 (6*) = constant H / a{r])d'r], (3-12) 

2 Jo 

Here, m is the classical mass of the particles involved in the scattering. Using the "semi- 
classical Levinson theorem" to determine the constant in this formula we deduce, using 
( pT9| ), the relation 

P+{e) = e'''^^^\ (3.13) 



where 

Jo tanhi(r7 + z^)tanhi(77-z^ 
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u I tanh^ r? tanh^ ^ 
2d{9)=2nB7r + - log ---——^J—-^^ (3.14) 



and ub is the number of bound states. Since according to the preceeding discussion 
nB = |: = p-forO<^< as /3 ^ we see that formula (|3.9|) is in complete 
agreement with the classical phase shift (|2.19|) and that ^ is proportional to (pQ. In the 
region — < ^ < there are no physical poles for P+ and ns = 0. A similar discussion 
can be carried out for ^ < and P_ . As /3 — > the number of physical poles of P_ varies as 
Ub = = and again we have agreement with the semi-classical Levinson theorem. 

From comparison of ( p.9|) and the classical phase shift we see that ^ and (fo are related 
linearly as ^ = ^fo- This leads correctly to the emission of a zero momentum soliton 
with the ground state degeneracy as discussed in 0. This linear relation can hold only 
for 1^1 < Beyond that value the ground state changes. To compare the quantum 

result with the classical phase shift we must then correlate appropriately the value of (p at 
infinity in the latter. The net effect is to replace (po by (po — 2tt. Eventually, the variation 
of ^ with (fo is therefore 

as illustrated in figure 4. It is very likely that (|3.15| ) is exactly true for finite f3 as well. 

We can finally recover ( |3.9| ) without appealing to our knowledge of quantum boundary 
bound states. For this one has to evaluate the action for the classical configuration. 
Following the discussion in @] we have, 

2m - C(,,H ^ I ^,log ^^^^ . ^ ^^^^ if 1/3.01 (3.16) 

where we used the soliton mass m = and 5 satisfies the differential equation 

P + OO / P+OO \ 

S{e) -tanhe6'{e) = J dti^j rfa;(^2 -Ssinh^tanh^j . (3.17) 

We restrict to the situation where we send in a kink, and (fQ is positive and so there is 
an anti-kink at the origin (see figure 1). In this case there are quantum boundary bound 
states. It is difficult to perform the double integral for the three-soliton solution explicitly 
because of the cumbersome expression for the integrand (Al). One might hope that in 
order to find the ^-independent piece of the phase shift it is sufficient to evaluate both 
sides of ( ^.17| ) in the limit as ^ — +oo or ^ — 0. However both of these limits do not 
seem to be helpful, because due to the non-uniform convergence of the right hand side 



(3.15) 
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of ( p.l7 ) it is not allowed to interchange such a limit with the integration. We evaluated 
the right hand side of ( |3.17] ) (with fixed 9) numerically for several different values of (j}Q 
using Mathematica. Combining ( p.l6| ) and ( |3.17|) , we obtain an estimate for C{(Pq), in 
agreement with C((/?o) = with accuracy 0.1%. This result was checked for different 

values of 9. We therefore obtain agreement with the semi-classical Levinson theorem. 

It is also interesting to investigate the cross-unitarity relation of in this semi- 
classical limit. In general we have 



~2 



ITT 



9 = St{29)P- [— + e]+ Sr{29)P+ {— + 9 



in 



(3.18) 



where S are the bulk S matrix elements. In the semi-classical limit we are considering, the 
reflection amplitude Sr vanishes. In this same limit we use ( p.9| ) to find 

P+{f-9) S{f-9) ^(^-9 + zi)s''- - 



P-{f + 9) S{f+9) 



ITT 

2 



9-i^ 



S(f + 9 + inS(^ + 9-i^ 



(3.19) 



' 5(i7r-26>) l'/' 
_S{i7r + 29)_ 

We simplify this expression by using the following relations deduced from unitarity and 
crossing symmetry at this order [|16| 



in 



St —-9 



S{9) 
S{9) 



S{-9) 

e ^"^5t(^), 



(3.20) 



to obtain ( p.l8|) indeed. 



4. Conclusions 

Although the method of images works nicely in the classical theory, it does not seem to 
extend to the quantum case: we have not been able to recast the boundary S matrix of 
as a product of bulk S-matrix elements. A related phenomenon is the non-trivial structure 
of the boundary S-matrix with one-loop corrections (the semi-classical expressions being 
the tree approximation). Recall that in the bulk one has 

5i(^;y) =exp j'dv In ts.nh' , (4.1) 
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where 



Scaling 

one finds the foUowing next-to-leading correction to the boundary S-matrix in the Dirichlet 
problem: 



^ , .„ ^ , „ „ tanh 1 — 

(4.3) 

In addition to the usual replacement of by we see the appearance of an entirely new 
factor involving the square-root of a hyperbolic tangent. 

It would be interesting to investigate further the structure of boundary bound states, 
most particularly in the quantum theory. Also, besides the semi-classical limit, one could 
also take the non-relativistic limit of the S-matrix. The appropriate candidate for the 
resulting non-relativistic potential would be probably the exactly-solvable Poschl- Teller 
model on half-line [p!7| : 

U[x) = Uo\——^ -2 , (4.4) 

V smh ax cosh ax / 

where a, /i, p and Uq are some constants, related to the parameters /3 and (po of sine- 
Gordon model. Note that (|4.4|) is simply a superposition of the repulsive hard-core kink- 
kink potential (first term) and the kink-antikink potential (second term), which are known 
from the bulk theory [^. This problem will be studied elsewhere [|T8| . 

Acknowledgements: 

We would like to thank P. Fendley for interesting discussions and A.Kapustin for the 
help with numerical computation and many useful remarks. N.W. would like to thank 
the Theory Division at CERN for its hospitality while this work was completed. This 
work was supported by the Packard foundation, the National Young Investigator program 
(NSF-PHY-9357207) and the DOE under grant number DE-FG03-84ER40168. 

Appendix 

We give here the explicit general form of classical solutions and directly extract the 
time delay (|2.19| ) in the limit M = oo. 
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In the center of mass reference frame the solution to ( |2.1|) reads: 



2e^'^^^~'^^ Ch(tsh9) ih e^""^ =P gx(l+2ch6') — 201—03+2 log(K ■;;) 
^ = ip4 arctan 2^x(l+che)-ai-a3+2log{u)^^fg-^()^ ^ g2xch6>-2ai+21og(?;) ;p ]^ ' 

where w = tanh(|), v = tanh(0), and we have set ai = 02, which means that the solution 
is invariant under the transformation t —t. The upper (resp. lower) sign corresponds 
to the situation when a stationary soliton (resp. anti-soliton) is used to adjust the value 
of field at the boundary. Solution (Al) refers to the case when the incoming and outgoing 
particle is the soliton with asymptotic value (p = 2it at x = +00. 

Let us represent (Al) in the form of rational function of variable ch(tsh^): 

(f)(x, t) = 4 arctan ^ — —r-} — r-Jr- 
' ' ri + sich(tsh^) 

The condition (j){x = 0, t) = ^0 implies that 

- = - = tanf^V {A2) 



ri si \ 4 

from which follows immediately 

1 /^^ 2, ^0 
03 = log I tan — 

For the argument of logarithm to be positive one should take ^0 > with the stationary 
anti-soliton and (j)Q < Q with the stationary soliton. This is illustrated in figure 1. Further, 
we obtain from (A2) 



2a 1 = log 



2 2 
u V 



w2 + tan2(^) 



(A3) 



1 + M2tan2(^) 

which agrees with ( p.l9|) where a = 2ai. Note that the time delay, obtained by (A3), is in 
fact always a time advance in both the attractive and repulsive cases. For the same value 
of ^0 the time delay for the soliton that lives on the left half-line a; < is not the same 
as that of the right half- line soliton (except for (^0 = ±7r). It differs by the sign of power 
factor in the formula ( |2.19| ) . The position of the "left" soliton is not the exact mirror image 
of the "right" soliton for generic (j)Q. 
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Figure captions 

Figure 1: a solution of the classical sine-Gordon equation with fixed boundary conditions. 
For (f)\x=o = 00 = ^ the solution is constructed out of a left-moving soliton, its right- 
moving image at a; < and the stationary anti-soliton in the middle (upper graph). For 
(l)\x=o = 00 = — ^ the solution is built out of three solitons (lower graph). 

Figure 2: a boundary bound state (boundary breather) for the values = ^ and 0o = x 
at four different instants of time. 

Figure 3: two possible ground state configurations for 00 = The configuration with 
asymptotic behaviour — > at infininty has lower energy than the other one. 

Figure 4: variation of ^ with 0o- 
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Figure 4: Variation of ^ function of (fo. 



